Analytic continuation of the critical line: suggestions for QCD 
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We perform a numerical study of tiie systematic effects involved in the determination of the critical 
line at real baryonic chemical potential by analytic continuation from results obtained at imaginary 
chemical potentials. We present results obtained in theories free of the sign problem, such as two- 
color QCD with finite baryonic density and three-color QCD with finite isospin chemical potential, 
and comment on general features which could be relevant also to the continuation of the critical line 
in real QCD at finite baryonic density. 
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I. INTRODUCTION 

The study of QCD at nonzero baryonic density by nu- 
merical simulations on a space-time lattice is plagued by 
the well-known sign problem: the fermionic determinant 
is complex and the Monte Carlo sampling becomes un- 
feasible. 

One of the possibilities to circumvent this problem is 
to perform Monte Carlo numerical simulations for imagi- 
nary values of the baryonic chemical potential, where the 
fermionic determinant is real and the sign problem is ab- 
sent, and to infer the behavior at real chemical potential 
by analytic continuation. 

The idea of formulating a theory at imaginary /i was 
first suggested in Ref. [ij, while the effectiveness of the 
method of analytic continuation was pushed forward in 
Ref. Since then the method has been extensively 
applied to QCD d i, i, i, 0, i, i, [l3l and tested in 
QCD-like theories free of the sign problem (Til . IT^ . ITsl . 
[l4 . [Tsl . [l6j and in spin models [13, [ll| • 

The state-of-the-art is the following: 

(i) the method is well founded and works fine within 
the limitations posed by the presence of nonana- 
lyticities and by the periodicity of the theory with 
imaginary chemical potential [19]; 

(ii) the analytic continuation of physical observables is 
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improved if ratios of polynomials (or Pade approx- 
imants [23|) are used as inter pola ting functions at 
imaginary chemical potential |l3l . [isj ; 

(iii) the analytic continuation of the (pseudo-)critical 
line on the temperature - chemical potential plane 
is well justified, but a careful test in two-color 
QCD |15|] has cast some doubts on its reliability. 

In particular, the numerical analysis in two-color QCD 
of Ref. [15] has shown that, while there is no doubt that 
an analytic function exists which interpolates numerical 
data for the pseudocritical couplings for both imaginary 
and real jj, across fi = 0, determining this function by an 
interpolation of data at imaginary fi could be misleading. 
Indeed, in the case of polynomial interpolations, there is 
a clear indication in two-color QCD that nonlinear terms 
in /i^ play a relevant role at real /i, but are less visible 
at imaginary /i, thus calling for an accurate knowledge of 
the critical line there and, consequently, for very precise 
numerical data. 

The above described scenario could well be peculiar 
to two-color QCD and strongly depends on the choice 
of parameters of Ref. p^ . Therefore, in this work we 
perform a systematic study of the analytic continuation 
of the critical line 

(i) in two-color QCD with a fermionic mass different 
than in Ref. [l5|, 

(ii) in another sign-free theory, SU(3) with a nonzero 
density of isospin. 

The aim of this study is to single out some general fea- 
tures of the analytic continuation of the critical line and 
to understand if and to what extent they can apply also 
to the physically relevant case of QCD. 
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FIG. 1: Susceptibility of Polyakov loop, chiral condensate, 
plaquette vs /3 in SU(2) on a 16^ x 4 lattice with am=0.2 and 
{afif = -0.09 (for the Polyakov loop), {afif = 0.0225 (for 
the chiral condensate) and (a/i)^ = —0.04 (for the plaquette). 
The solid lines represent the Lorentzian interpolation. 



II. ANALYTIC CONTINUATION OF THE 
CRITICAL LINE IN TWO-COLOR QCD 

As in Ref. [13], to which we refer for further de- 
tails about our numerical setup, we have performed 
Monte Carlo simulations of the SU(2) gauge theory with 
Nf = 8 degenerate staggered fermions. The algorithm 
adopted has been the usual exact algorithm described 
in Rcf. 21], properly modified for the inclusion of a finite 
chemical potential. 

The general strategy is the following: first the critical 
line is sampled by determining, for a set of values, the 
critical couplings /3c (z^^) for which the susceptibilities of 
a bulk observable exhibit a peak. Since the location of 
the critical line may well depend on the chosen observ- 
able, we have used three probe observables: the chiral 
condensate, the Polyakov loop and the plaquette. The 
/i^ < values are chosen to lie inside the so-called first 
Roberge- Weiss (RW) sector, i.e. in the range delimited 
by (afi)^ = -(7r/8)2 and /i^ = (see Ref. 13] for a 
detailed discussion on the phase diagram in the T — ifi 
plane) . 

Then, the critical line is guessed by interpolating the 
values of PdfJ-'^) for /i^ < only. The validity of the inter- 
polation is evaluated by comparing its analytic continua- 
tion to the region /x^ > with the direct determinations 
of the critical coupling in this region. 

The procedure above is repeated for each of the three 
probe observables considered. 



A. Nt = 4, am = 0.07 

The SU(2) gauge theory with Nf — 8 degenerate stag- 
gered fermions of mass am = 0.07 has been studied on a 
16^ X 4 lattice in Ref. [H]. We found that 

(i) there is a very weak dependence of the values of 
/3c(m^) on the probe observable; 

(ii) there is no better interpolation of /3c (a*^) data at 
/i^ < than a polynomial of the form ag + ai(a/i)^; 

(iii) the extrapolation to the region jj? > overshoots 
the direct determinations of /3c (m^) i'^ that region, 
the discrepancy becoming larger and larger for in- 
creasing n"^; 

(iv) data for /3c (^^) for both /i^ < and /i^ > can 
be nicely interpolated by an analytic function (a 
polynomial of third order in works nicely), this 
excludes the option that the critical line be not 
smooth and possibly nonanalytic in the thermody- 
namic limit. 

The conclusion we drew in Ref. [l^ was that terms 
of order /x^ and play a relevant role at /i^ > 0, but 
are less visible at /x^ < 0, thus calling for an accurate 
knowledge of the critical line in all the first RW sector. 
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TABLE I: Summary of the values of fic{p?) obtained by fitting the peaks of the susceptibilities of chiral condensate {ipip), 
Polyakov loop (L) and plaquette (P) in SU(2) on a 16"^ x 4 lattice with fermionic mass am—Q.2. For each interpolation the 
X^/d.o.f. is given. 





{^^) 


x7d.o.f. 


{L) 


xVd.o.f. 


{P) 


xVd.o.f. 


-0.1452 


1.8098(71) 


1.01 


1.7991(14) 


0.51 


1.8058(67) 


1.13 


-0.1225 


1.7713(30) 


0.97 


1.7709(30) 


0.21 


1.7688(43) 


0.20 


-0.09 


1.7465(37) 


0.30 


1.7365(76) 


0.66 


1.7448(60) 


0.79 


-0.04 


1.6998(42) 


0.92 


1.6934(53) 


2.10 


1.6948(86) 


0.73 


-0.01 


1.6762(41) 


0.37 


1.6680(86) 


0.20 


1.6745(77) 


0.24 


0. 


1.6694(44) 


0.30 


1.6634(76) 


1.50 


1.6709(71) 


0.58 


0.0225 


1.6445(20) 


0.03 


1.6180(65) 


1.01 


1.6440(34) 


0.14 


0.04 


1.6430(54) 


1.04 


1.6105(48) 


1.87 


1.6367(67) 


0.74 


0.09 


1.5989(45) 


0.30 


1.5606(70) 


0.73 


1.6005(81) 


0.68 



TABLE II: Parameters of the fits to the critical couplings 
in SU(2) on a 16^ x 4 lattice with fermionic mass am=0.2 
according to the fit function /3c (/x^) = ao + ai(o/x)^, for each of 
the three probe observables used (chiral condensate, Polyakov 
loop and plaquette). 



observable 


ao 


ai 


xVd.o.f. 


(^^> 


1.6671(28) 


-0.876(32) 


1.61 


{L) 


1.6556(44) 


-0.981(33) 


1.34 


(P) 


1.6656(48) 


-0.885(49) 


1.40 



Verifying if this scenario is peculiar to SU(2) with the 
parameters choice as in Ref. UM is the main motivation 
of the present work. For this reason, in the following 
subsection we repeat the same analysis of Ref. [ij] in 
SU(2) with a different numerical set up, whereas in the 
next section we turn to a completely different theory. 



B. Nt ^4, am = 0.2 

We have considered SU(2) with Nf = 8 degenerate 
staggered fermions of mass am — 0.2 on a 16'^ x 4 lattice. 
The reason for setting a larger value for the quark mass 
is simple: at any fixed temperature T the critical value of 
fi gets larger, so that the critical line could become less 
curved in the physical region /i^ > 0. This should reduce 
the importance of higher orders in /i^ in the description 
of the critical line by a polynomial. 

As in Ref. the critical line is sampled by looking 
for peaks in the susceptibilities of Polyakov loop, chiral 
condensate and plaquette. In Fig. [1] we show as exam- 
ple the behavior in f3 of the susceptibility of the three 
observables at different values of /x^, together with the 
Lorentzian interpolation of the peak. The statistics for 
each data point is about 10-20k trajectories of one Molec- 
ular Dynamics (MD) time length. 

All determinations of /3c (/i^) are summarized in Ta- 
ble [H from which we see a mild dependence on the probe 
observable for /z^ < 0; at /x^ > 0, instead, determinations 
from the Polyakov loop are systematically below those 
from the chiral condensate and the plaquette, which in- 



stead agree within errors. 

As in Ref. [l^, the best interpolation of /3c (m^) data 
at /i^ < is a polynomial of the form op -I- ai(a/i)^. In 
Table |lT] we report the values of the fit parameters and 
their uncertainties as obtained with the MINUIT min- 
imization code. The extrapolation to /x^ > compares 
very well with the direct determinations of /3c (/i^) in that 
region, as shown in Figs. [H In these figures, and in the 
following of the same kind, the band around the best fit 
(dotted line) represents the 95% confidence level (CL) 
region. 

We can draw some partial conclusions about the ana- 
lytic continuation of the critical line in SU(2): 

• for our given accuracy, there seems to be no room 
for quartic and sextic terms in /x^ in the interpola- 
tion of /3c(/x^) data for imaginary /i; 

• the extrapolation to /i^ > works definitely better 
for larger masses, i.e. away from the chiral limit. 

Both these features could be peculiar to the SU(2) the- 
ory. This motivates the need to repeat a similar analysis 
in a different theory, sharing with SU(2) the property of 
being free of the sign problem. 

III. ANALYTIC CONTINUATION OF THE 
CRITICAL LINE IN THREE-COLOR QCD AT 
FINITE ISOSPIN CHEMICAL POTENTIAL 

The SU(3) gauge theory with a finite density of isospin 
charge [1^ is a theory in which the chemical potential is 
fiiso for half of the fermionic species and — /iiso for the 
other half. The partition function, which is even in /iiso 
and depends only on /ifgQ, can be written as follows: 

Z{T, /iiso) = J VUe~^^ det M[^l,,,] det M[-^l,,,] (1) 

where the integration is over gauge link variables, Sg is 
the pure gauge action and M the fermion matrix (we 
adopt a standard staggered discretization). This leads 
to a real and positive measure, because of the property 
det Af[— /iiso] = (det M [/iiso])*, and therefore to a theory 
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FIG. 2: Critical couplings obtained from the susceptibilities 
of Polyakov loop, chiral condensate and plaquette in SU(2) 
on a 16^ X 4 lattice with am=0.2, together with a linear fit 
(dotted line) in (a/x)^ to data with /x^ < 0. Here and in the 
following plots of the same kind, the solid lines around the 
best fit line delimit the 95% CL region. 



free of the sign problem. This theory is obviously closer 
to real QCD than two-color QCD, being yet unphysical, 
since it implies a zero baryonic density, while in nature 
a nonzero isospin density is always accompanied by a 
nonzero baryonic density; moreover the isospin charge is 
not a conserved number in the real world. 

Nevertheless, for our purposes this theory is very con- 
venient since it provides us with another theoretical lab- 
oratory for the method of analytic continuation. 

Similarly to SU(2) with finite baryonic density, at 
imaginary values of the chemical potential fiiso (which 
we will call simply from now on) the theory exhibits 
RW-like transition lines (we refer to the following sub- 
section for a discussion about the determination of their 
position). The first RW sector is given by the strip 
-(0.5)2 < ^y(T,T)^ < 0. 

In our numerical analysis, we consider finite isospin 
SU(3) with Nf — 8 degenerate staggered fermions of 
mass am = 0.1 on a S'^ x 4 lattice. Differently from SU(2) 
with finite baryonic density, the critical line in the tem- 
perature - chemical potential plane is a line of (strong) 
first order transitions, over all the investigated range of 
/i^ values, -0.2304 < n'^/inTf < 0.2025. This is one 
of the reasons for working on a smaller volume than in 
the SU(2) case (tunneling between the different phases 
would have been sampled with much more difhculty on 
a larger volume) and clearly emerges from the distribu- 
tion on the thermal equilibrium ensemble of the values of 
observables like the (real part of) the Polyakov loop, the 
chiral condensate, and the plaquette across the transition 
(see Fig. [3] as an example showing the typical two-peak 
structures) . As a further evidence of the first order na- 
ture of the transition, we show in Fig.[4]the Monte Carlo 
run history of the (real part of) the Polyakov loop at 
^i'^|{■KTf = -0.09 and (3 = I3c = 4.750, which exhibits 
tunneling events between the two phases every few thou- 
sands trajectories, on average. Typical statistics have 
been around lOK trajectories of 1 MD unit for each run, 
growing up to lOOK trajectories for 2-3 /3 values around 
/3c (m^), for each /x^, in order to correctly sample the criti- 
cal behavior at the transition. The critical Pifi^) is deter- 
mined as the point where the two peaks have equal height 
and, in all the cases considered, this point turned out to 
be the same for all the adopted observables. In particular 
this means that, at the transition, both the chiral con- 
densate and the (real part of) the Polyakov loop "jump" 
at the same point (see Fig.[Sl as an example). 

If one looks at the "gap" between the two peaks at 
the transition (determined as the difference between the 
values of an observable at the maxima of the right and 
left distribution peaks), it turns out that the depen- 
dence of this gap on /^^/(TrT)^ in the considered range 
is slightly observable-dependent, being roughly constant 
for the Polyakov loop and increasing for the chiral con- 
densate and the spatial plaquette, up to a turning point 
at jj? /{ttTY — 0.1, where the strength of the transition 
seems to reach a maximum. In particular, no critical 
endpoint for the first order line is detectable, within the 
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FIG. 3: Distributions of the real part of the Polyakov loop (left) and of the chiral condensate (right) in SU(3) with finite isospin 
density on a 8'^ x 4 lattice with am— 0.1 at fi'^ /(ttT)^ =0.16 and for three fJ values around the transition. 
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FIG. 4; History of the Monte Carlo run of the real part of the 
Polyakov loop in SU(3) with finite isospin density on a 8^ x 4 
lattice with am=0.1 at fi'^/iirTf = -0.09 and /3=4.750. 



FIG. 5: Real part of the Polyakov loop and chiral condensate 
vs f3 in SU(3) with finite isospin density at fj? /(iiT)^ = —0.09 
on a 8^ X 4 lattice with am=0.1. 



A. RW-like transitions at finite isospin density 



explored range. We notice that in Ref. [1^, where the 
same theory has been explored even if with a slightly 
larger quark mass {am — 0.14), a critical endpoint was 
reported, but for a real chemical isospin chemical poten- 
tial {n/T ^ 2.5) which is well outside the range explored 
by us. 

When combining the first RW-like transition line 
with the numerical determinations for the critical /3's 
(see Sec. IIII B[) . the resulting phase diagram in the 
(/3,AtV(7rr)2) looks like in Fig. [3 



In order to understand the unphysical, RW-like phase 
transition taking place in the high temperature region 
of the QCD phase diagram in the presence of an imag- 
inary isospin chemical potential, let us first remember 
the origin of the usual RW transition in the presence of 
an imaginary baryon chemical potential. That is linked 
to the dynamics of the Polyakov loop. Let us define 
9 = Im(/z)/T, where fj, is the baryon chemical potential, 
and let us call e'^J , with j — l,...Nc, the eigenvalues 
of the untraced Polyakov loop, which satisfy the condi- 
tion '^j(t>j = Oinod27r- The onc-loop effective potential 
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FIG. 6: Gap of several observables at the transition for several 
values of the isospin chemical potential. 
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FIG. 7: Sketch of the phase diagram of SU(3) with finite 
isospin density on a 8^ x 4 lattice with am=0.1. 



of the untraced Polyakov loop in presence of the baryon 
chemical potential is then given by [l^ 



for = Vj, and the center element corresponding to a 
real Polyakov loop is selected. For 9^0 the minimum of 
the fermion contribution moves and as 9 crosses t:/Nc the 
absolute minimum of the potential moves to a different 
center element (RW transition). 

In the case of an isospin chemical potential instead, 
the chemical potential is fi for half of the fcrmionic fla- 
vors and —fi for the other half. The effective potential 
therefore looks as follows: 




( 




0/c 


-0 






TT 


mod2 / 



mod2 



- 1 



3(12 



(3) 



In this case the fermion contribution is left unchanged 
as — > —(pj, i.e. it is symmetric under complex conju- 
gation of the Polyakov loop. However it easily verified 
that as 9 crosses a critical value 9c — 0.5168 tt, the two 
degenerate minima corresponding to complex center el- 
ements becomes deeper than then one corresponding to 
a real center element, hence for 9 > 9c the conjugation 
symmetry of the Polyakov loop is spontaneously broken 
and one of the two complex minima is selected. 

Notice that the above effective potential is computed 
at one loop and for zero mass fermions. Higher order and 
finite quark mass corrections do not affect the location 
of the RW transition in the case of an imaginary baryon 
chemical potential, which is fixed at 6* = {2k + l)Tr / Nc by 
symmetry. The same corrections may however influence 
the location of 9c for isospin chemical potential, hence 
9c may be quark mass and temperature dependent. In- 
deed, as shown in Fig. [7l numerical simulations locate 
approximately 9c ~ 0.48 n as the point where the un- 
physical RW-like line meets the analytic continuation of 
the physical transition line. 

A more detailed analysis of the QCD phase diagram in 
the presence of imaginary baryon and isospin chemical 
potentials goes beyond our present purposes and will be 
presented elsewhere [23 ]. 





(f>j <t>k 



inod2 




The first term is the gluonic contribution, which has 
three degenerate minima, corresponding to SU(iVc) ma- 
trices belonging to the center of the group. In absence of 
fermions that means spontaneous breaking of center sym- 
metry. For 9 = Q the fermion contribution is minimum 



B. Results for the critical line at finite isospin 

In Table IIIII we summarize our determinations of the 
critical couplings for the finite isospin SU(3) theory on a 
8'^ X 4 lattice with fermionic mass am^O.l. 

We observe from the very beginning that data for 
/3c (m^) for both < and /i^ > can be globally fitted 
by an analytic function (a polynomial of third order in 
^ jiirTY nicely works). The fit parameters of the global 
fit are given in the last line of Table IIVI while Fig. [5] 
shows how the fit compares with the data. 
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FIG. 8: Critical couplings obtained in SU(3) with finite 
isospin density on a 8'' x 4 lattice with am=0.1, together with 
a polynomial fit of order /i® to all data. 



The question is if there are interpolations of the critical 
couplings at /i^ < only, that, when continued to > 0, 
agree with the critical couplings directly determined in 
the latter region. 

We have tried several kind of interpolations of the crit- 
ical couplings at /i^ < 0. At first, we have considered 
interpolations with polynomials up to order /i^" (see Ta- 
ble [IV] for a summary of the resulting fit parameters). 
We can see that data at < are precise enough to be 
sensitive to terms beyond the order indeed, a good 
X^/d.o.f. is not achieved before including terms up to the 
order /i^, in agreement with the outcome of the global fit 
discussed above. The extrapolation to /z^ > for the 
polynomial of order /i^ is shown in Fig. [5] (right), in com- 



parison with that obtained from the polynomial of order 
fi^ (Fig. [H] (left)). The extrapolation agrees with direct 
determinations of Pcip"^), within the 95% CL band, only 
if nonlinear terms in /i^, at least up to /i^, are taken into 
account. 

Then, we have considered interpolations with ratios of 
polynomials of order up to /i® (see Table |V] for a sum- 
mary of the resulting fit parameters). In all but one case 
we got good fits to the data at /i^ < 0, but only two ex- 
trapolations to /i^ > compare well with numerical data 
in that region: the ratio of a 4th to 6th order polyno- 
mial and the ratio of a 6th to 4th order polynomial (see 
Fig. [TU|) . It is interesting to observe that the two inter- 
polations which "work" have in common the number of 
parameters. 

Both kinds of fits considered so far have evidenced that 
the role of terms of order larger than cannot be ne- 
glected. Since the data more sensitive to these terms are 
those farther from fj? = 0, while data closer to — 
should be enough to fix the coefficient of the fi^ term in 
a polynomial interpolation, we followed a different strat- 
egy: we performed a fit with a polynomial of the form 
ao + aifiyinT)^ in the region nU{nT)^ < iJ^/{T:Tf < 
and varied /^^j,; in order to find the largest possible 
interval in which the fit gives both a good x'^/d.o.f. and 
a stable value (within errors) for the linear coefficient ai; 
this fixes the value of the parameter ai with some accu- 
racy (see Table IVl)) ^. After that, we fixed ai at 0.590 
(allowing an arbitrary variability of 0.006) and repeated 
the fit on all available data at < with a polynomial 
of the form uq + aifj."^ /{nTy + a2At^/(7rT)4 + a3fj.^/{nT)^. 
The resulting interpolation and its extrapolation to the 
region /i" > are shown in Fig. 1111 The comparison with 
critical couplings at //-^ > is good and the 95% CL band 
is narrower than in the unconstrained /i^-polynomial fit 
(see Fig. ini (right)), meaning that this procedure leads to 
increased predictivity for the method of analytic contin- 
uation. 

At last, we have attempted a fit strategy never used 
so far in our studies on the present subject and, to our 
knowledge, nowhere in the literature on the analytic con- 
tinuation from imaginary chemical potential. The idea 
is to write the interpolating function in physical units 
and to deduce from it the functional dependence of Pc 
on /i^, after establishing a suitable correspondence be- 
tween physical and lattice units. The natural, dimen- 
sionless variables of our theory are T/Tc{0), where Tc(0) 
is the critical temperature at zero chemical potential, and 
fi/^TrT). The question that we want to answer is if fit- 
ting directly the dependence of T/Tc{0) on ijl/{ttT) may 
lead to increased predictivity for analytic continuation. 
We shall name this kind of fits, made directly in terms 
of the dimensionless physical variables of the theory, as 
"physical" fits. 



^ A similar strategy has been followed also in Ref. [25 
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TABLE IV: Parameters of the fits to tlie critical couplings in SU(3) with finite isospin density on the 8'^ x 4 lattice with fermionic 
mass am=0.f according to the fit function /3c(At^) = ao + ai^^ /{ixTf +a2^i'^ /{i^Tf + a3^i^ /{iiTf + aA^l^ /{-ivTf + as^^° /{iiTf^ . 
Blank columns stand for terms not included in the fit. The last column reports the largest value of /{■rTY' included in each 
fit. The * in the last but one line means "constrained parameter". 





ai 


02 


a-j, 


04 


as xVd.o.f. Mmax/(7rT)^ 


4.69535(79) 


-0.6153(53) 








6.33 


0. 


4.6982(13) 


-0.529(26) 


0.37(11) 






2.24 


0. 


4.6967(17) 


-0.636(72) 


-0.83(74) 


-3.4(2.1) 




0.88 


0. 


4.6966(19) 


-0.638(63) 


-0.81(32) 


-3.0(1.2) 


1.3(5.0) 


0.97 


0. 


4.6968(17) 


-0.621(24) 


-0.50(26) 


-1.54(45) 


1.00(53) 


-10.(19.) 0.94 


0. 


4.6974(8) 


-0.596(12)* 


-0.43(14) 


-2.36(61) 




1.07 


0. 


4.6970(10) 


-0.585(12) 


-0.190(44) 


-1.528(35) 




1.10 


0.2025 




FIG. 9: Critical couplings obtained in SU(3) with finite isospin density on a 8'^ x 4 lattice with am=0.1, together with a 
polynomial fit of order /i^ (left) and (right) to data with < 0. 



While fj,/{TTT) is one of the dimensionless variables 
used in our simulations, T/Tc{0) is not and must be de- 
duced from the relation T = l/{Nta{(3)), where Nt is the 
number of lattice sites in the temporal direction and a(/3) 
is the lattice spacing at a given (3 ^. Since our determi- 
nations for /3c range between ~ 4.5585 and ~ 4.842 (see 
Table IIIip it can make sense to use for a(/3) the pertur- 
bative 2-loop expression with TVc = 3 and Nf = 8. 

We have tried several different fitting functions and 
report two cases which work particularly well. The first 
is given by the following 3-parameter function: 



Tc{0) 



1 + B^xy{7:TM)^ 
1 + A^l^/{^IT,{^l)Y 



(4) 



leading to the following implicit relation between 13c and 



^ Strictly speaking the lattice spacing depends also on the bare 
quark mass, which in our runs slightly changes as we change /3 
since we fix am. However in the following evaluation, which is 
only based on the perturbative 2-loop /3-function, we shall neglect 
such dependence. 



2. 



a(/9c(M'))'l 2-loop = a(/9c(0))'|2-ioop 



(5) 



In 



Fig. 

,2\ 



[T^ left) we compare 



the values of 

a(/3c(Ai'^))'^ I |2-ioop obtained using our determina- 
tions for Pdn^) with the interpolation according to the 
fit function ([S]) to data in the region /i^ < 0: one can see 
that the extrapolation to the region /i^ > behaves very 
well. In Figs, [l^right') we have rescaled the vertical axis 
in order to give the behavior of /5c (/^^) versus the isospin 
chemical potential. The values of the fit parameters are 



/?c(0) = 4.6977(13) , 



A = -3.25(26) , 
B = -2.62(12) , 



(6) 



with x^/d.o.f.=1.33. We observe that /3c(0) is in nice 
agreement with the direct determination (see Table |TTT]). 

As an alternative function for the shape of the critical 
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TABLE V: Parameters of the fits to the critical couplings in SU(3) with finite isospin density on the 8^ x 4 lattice with 
fermionic mass am=0.1 according to the fit function f3c{tJ,^) = (ao + ai/i^/(7rr)^ + 02^'*/ {■kT)'^ + a3fi' /{TvT)^)/{l + a4fj,^ / (vrT)^ + 
aryfi'^ / (tvT)* + aefj,^ /(ttT)^). Blank columns stand for terms not included in the fit. 



ao 


0.1 


0,2 


as 


04 


as 


ag xVd-o.f. 


4.6982(13) 


2.21(76) 






0.58(16) 




2.04 


4.69748(79) 


18.3737(83) 






4.0328(17) 


0.4711(58) 


0.50 


4.6970(12) 


15.034(17) 






3.3287(35) 


0.461(17) 


0.236(55) 0.49 


4.69730(85) 


17.3130(99) 


-2.177(33) 




3.8089(21) 




0.44 


4.6967(14) 


2.124(34) 


-20.38(17) 




0.5855(73) 


-4.136(35) 


0.74 


4.6971(13) 


1.920(18) 


-58.143(93) 




0.5353(38) 


-12.277(19) 


-1.332(67) 0.55 


4.6971(11) 


16.284(13) 


-2.202(62) 


-0.46(20) 


3.5928(27) 




0.46 


4.6971(13) 


4.503(19) 


-44.770(95) 


4.82(33) 


1.0851(39) 


-9.366(20) 


0.54 


4.6971(14) 


4.956(19) 


-42.960(97) 


4.55(40) 


1.1816(39) 


-8.969(20) 


-0.007(82) 0.65 




FIG. 10: Same as Fig. |9l with the fitting function given by the ratio of a 4th to 6th order polynomial (left) and the ratio of a 
6th to 4th order polynomial (right). 



4.9 



4.8 



4.6 



4.5 





1 1 1 1 1 


1 1 1 1 1 

SU(3) finite isospin 












polynomial order jj.* 






constrained 




I.I. 


1 . iV\ 



0.1 



|X^/(7UT)^ 



FIG. 11; Same as Fig. [9] with the fitting function given by a 
polynomial of order fJ" with the quadratic term constrained 
as explained in the text. 



TABLE VI: Values of the parameter ai resulting from the fit 
to the critical couplings in SU(3) with finite isospin density 
on the 8'^ X 6 lattice with fermionic mass am=0.1 according 
to the fit function /3c(a'^) ~ ao + aifi^ /{tvT)'^ , in the interval 

lfiLn/{^Tf,0]. 





ai 


xVd.o.f. 


-0.04 


-0.60(5) 




-0.0625 


-0.58(3) 


0.39 


-0.09 


-0.585(21) 


0.29 


-0.126025 


-0.586(12) 


0.22 


-0.16 


-0.5915(88) 


0.39 


-0.1849 


-0.5931(77) 


0.39 


-0.198025 


-0.5945(69) 


0.39 


-0.2116 


-0.6010(62) 


1.83 


-0.2209 


-0.6078(57) 


3.72 



line, we have also tried the following 




-A) [cos(^)]'' 
{I- A) [cosh(f)]- 



(7) 



10 




-0.2 -0.1 0.1 0.2 

^'/(TlT)' 



4.9 



4.8 



P 4.7 



4.6 



4.5 








1 1 1 1 1 
SU(3) finite isospin 






8^x4 


- 




"physical" fit 




1,1, 


1 , 1 \,\\ 



-0.2 -0.1 0.1 



0.2 



FIG. 12: Values of [a{l3c(}i^))kf (left) and (right) in SU(3) with finite isospin density on a 8^ x 4 lattice with am=0.1, 

together with the fit to data with /i^ < according to the fit function ([5]). 
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"physical" fit 



_L 




_L 



-0.2 -0.1 0.1 0.2 

lx'/(7iT)' 



4.9 



4.8 



P 4.: 



4.6 



4.5 





1 1 1 1 1 


1 1 1 1 1 
SU(3) finite isospin 






8^x4 






"physical" fit 




1,1, 


1 , 1 \-\ 



-0.1 0.1 



FIG. 13: Same as Fig. [12] with the fit to data with <Q according to the fit function ((T]). 



which expHcitly encodes the expected periodicity of the 
partition function for imaginary fj,. The fit to data at 
imaginary /i is very good and its extrapolation to the 
real chemical potential side compares impressively well 
with data (see Fig. [T3|) . The resulting fit parameters are 

/3^(0) = 4.6969(12) , A ^ 1.508(15) , 

B = 0.560(32) , (8) 

with x^/d.o.f.=0.39. This function is a good candidate 
to parametrize the critical line for small values of /i/T. 

In both cases, Eqs. ^ and [3 the physical fit has 
worked very well and with a reduced number of parame- 
ters with respect to our previous fits, leading to increased 
predictivity and consistency with data at real chemical 
potentials. In both cases one can easily check that the 



adopted functions are not appropriate for a continuation 
of the critical line down to the T = axis, but this is not 
the aim of our study since such extrapolation would be 
questionable anyway. 

We have used other functional forms for Tc{^), inspired 
by the Gross-Neveu model [2^ and by the random matrix 
theory [131 , and found that they allow one to interpolate 
data at /x^ < 0, but their extrapolation does not match 
data at /i^ > 0. 

Finally, in order to compare the efficiency of the dif- 
ferent predictions of the critical line, we present in Ta- 
ble lVIII the extrapolated value of /3c at a common value of 
the real chemical potential (/i^/(7rT)^ = 0.2) for all the 
interpolation methods adopted. We can see that the last 
two methods listed (constrained and physical fit) have 
a narrower error band, the physical fit having a better 
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agreement with data. 

TABLE VII: Extrapolated values of at ^?/{-KTf = 0.2 
for the interpolation methods adopted in Figs. [9lll3l 



/3c kind of interpolation 



4.51(6) 


Fig 


. IHtright) 


4.55(3) 


Fig 


. [inileft) 


4 'Sfi+° °^ 

^•3"-0.02 


Fig 


.Upright) 


4 54+0 02 

^■^4-0.01 


Fig 


m 


^■3' -0.02 


Fig 


. [m (right) 


4.56(1) 


Fig 


. Ha (right) 



IV. DISCUSSION 

The present study is part of a larger project which, 
based on the investigation of QCD-like theories which are 
free of the sign problem, is aimed at testing the validity 
of the method of analytic continuation and at improving 
its predictivity, in view of its application to real QCD. 

In particular, we have presented results concerning the 
analytic continuation of the critical line in 2-color QCD 
and in QCD with a finite density of isospin charge. We 
have detected some common features and developed some 
strategies, which could apply and be useful also for real 
QCD at finite baryon density. Let us briefly summarize 
them. 

(i) Non-linear terms in the dependence of the pseu- 
docritical coupling (3c on /i^ in general cannot be 
neglected. A polynomial of order /i^ seems to be 
sufficient in all explored cases. The prediction for 
the pseudocritical couplings at real chemical poten- 
tials may be wrong if data at imaginary /x are fitted 
according to a linear dependence. 

(ii) The coefhcients of the linear and non-linear terms 
in /i^ in a Taylor expansion of /3c (/i^) are all neg- 
ative. That often implies subtle cancellations of 
nonlinear terms at imaginary chemical potentials 
{jj? < 0) in the region available for analytic con- 
tinuation (first RW sector). The detection of such 
terms, from simulations at /i^ < only, may be 
difficult and requires an extremely high accuracy. 
As a matter of fact, the simple use of a sixth order 



polynomial to fit data at imaginary fi leads to poor 
predictivity, which is slightly improved if ratio of 
polynomials are used instead. 

(ii) An increased predictivity is achieved if the follow- 
ing strategy is adopted: first, the linear term in fj,^ 
is fixed from data at small values of /i^ only. Then 
a nonlinear fit (e.g. with a sixth order polynomial) 
is performed, keeping the term constrained. 

(iii) We have proposed a new, alternative ansatz to 
parametrize the critical line directly in physical 
units in the T, fj, plane (instead than in the /3, fi 
plane) and given two explicit realizations. We have 
shown that this "physical" ansatz provides a very 
good description of the critical line, moreover with 
a reduced number of parameters, and leads to an 
increased predictivity, comparable to that achieved 
by the "constrained" fit. 

We plan to apply our experience to the determination 
of the critical line for real QCD in the near future. The 
lesson from the present study is that the preferred extrap- 
olation methods to be considered in that case are a non- 
linear polynomial, at least of sixth order in /i, in which 
the /i^ term is constrained, or rather the new "physical" 
ansatz for the critical line that we have proposed. The 
two methods have shown comparable predictivity. The 
former requires dedicated simulations at small chemical 
potentials, in order to better fix the fj,"^ term, which could 
also be fixed by using alternative methods like reweight- 
ing or Taylor expansion. The latter is more appealing 
since it describes the critical line with a reduced number 
of parameters: this is very likely at the origin of its in- 
creased predictivity and is probably related to the fact 
that one is looking for a direct relation between two phys- 
ical quantities (T/Tc(0) and /i/T) in this case, without 
going through bare lattice quantities. 
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